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Abstract: We analyze the alignment of the low multipoles (quadrupole and octupole) of 
various maps of the WMAP 5yr release: the CMB maps obtained with ILC and MCMC 
methods, the CMB map in the V band after foreground reduction, and, for comparison, 
the (not cleaned) V band map. We study how much this alignment is polluted by residuals 
on the Galactic region. Among the considered maps, the WMAP-ILC turns out to be the 
most clean map from the point of view of the proposed test. This result has been found 
studying the redistribution (due to the masking process) of each bin of the probability 
distribution functions of the alignment estimators. By construction, our method, feasible 
through Monte Carlo simulations, works for any possible mask adopted in the analysis of 
data from current and forthcoming CMB anisotropy experiments and it can only exclude 
that the considered map is clean. 
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1. Introduction 

The anisotropy pattern of the cosmic microwave background (CMB), obtained by Wilkinson 
Microwave Anisotropy Probe (WMAP), probes cosmological models with unprecedented 
precision (see [1, 2] and references therein). Although WMAP data are largely consistent 
with the concordance A cold dark matter (ACDM) model, there are some interesting de- 
viations from it, in particular on the largest angular scales. They can be divided in the 
following categories. 1 ) Lack of power at large scales. The angular correlation function is 
found to be uncorreleted (i.e. consistent with 0) for angles larger than 60 degrees. In [3, 4] 
it has been shown that the probability associated to this event is low as 0.15%. Still in this 
category we mention the surprisingly low amplitude of the quadrupole term of the angular 
power spectrum (APS), already found by Cosmic Background Explorer (COBE) [5, 6], and 
now confirmed by WMAP [1, 2]. 2) Unlikely alignments of low multipoles. An unlikely 
(for a statistically isotropic random field) alignment of the quadrupole and the octupole 
is described in reference [7, 8, 9, 10, 11, 12]. Moreover, both quadrupole and octupole 
align with the CMB dipole [4]. Other unlikely alignments are described in [13, 14, 15]. 

3) Hemispherical asymmetries. It is found that the power coming separately from the two 
hemispheres (defined by the ecliptic plane) is too asymmetric (especially at low i) [16, 17]. 

4) Cold Spot. In Ref. [18] it is reported a detection of a non Gaussian behaviour in the 
southern hemisphere with a wavelet analysis technique (see also [19]). 

It is still unknown whether these anomalies come from fundamental physics or whether 
they are the residual of some not perfectly removed astrophysical foreground or systematic 
effect [20, 21, 22, 23, 24]. As an example of the latter kind, in references [25, 26] it is 
presented a study about the impact of the dipole straylight contamination on the low 
amplitude of the quadrupole and on the low £ alignments for Planck ^ characteristics and 
^http://www.rssd.esa.int/planck 
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capabilities. Many efforts liave been dedicated to the development of methods aimed at 
the discrimination between spurious and cosmological effects [27, 28, 29, 30]. The still 
open question about the origin of such anomalies has attracted a lot of interest in the 
last few years. In [31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41] there are some papers about 
the possible generation of the low ^ anomalies in the context of some specific models 
of fundamental physics. In [42] it is claimed that no model of contamination that is 
statistically independent of the source of the primary CMB anisotropy, can explain this 
large-scale power deficit. In other words, if a contamination (not taken into account) 
is responsible for the lack of power it must have a correlation with the primary CMB 
anisotropy. 

In the present paper we focus on the second of the aforementioned list of anomalies, 
i.e. on the unlikely alignments of low multipoles (quadrupole and octupole). We take 
into account various maps of the WMAP 5 year release: the ILC map, the MCMC map, 
the (not cleaned) V band map and the foreground reduced V band maps (see Fig. 1). 
The ILC (Internal Linear Combination) map, available at pixel size of 6.87 arcmin, has 
been obtained from a weighted linear combination of the five intensity maps at the various 
WMAP frequencies smoothed to 1 degree (FWHM) resolution, with weights chosen to 
maintain the CMB anisotropy signal while minimizing the Galactic foreground contribution 
in different 12 regions covering the whole sky. The WMAP team believes that it is suitable 
for analyses on angular scales greater than about 10 degrees. The MCMC (Monte Carlo 
Markov Chain) map is again smoothed to 1 degree resolution but with a pixel size of 
54.97 arcmin. Prom the MCMC fit the WMAP team derived anisotropy maps for the 
Galactic foreground diffuse components and that for the CMB component, considered in 
this work. We exploit also the foreground reduced (i.e. the CMB anisotropy map derived 
with a subtraction of the foreground components using a Foreground Template Model) 
map in the V band, the WMAP frequency channel where the original level of foreground 
is minimum. For comparison, we consider also the map in the V band without any kind 
of foreground subtraction. These two last maps centred at 61 GHz have the original beam 
resolution of 19.8 arcmin and are provided at a pixel size of 6.87 arcmin. Clearly, since 
we will analyse the large scale properties of the sky, the fact that the above maps have 
been provided at different resolutions and pixel sizes do not have any significant impact for 
this work. All these products are publicly available at LAMBDA web site ^ where further 
information can be found. 

In this paper we address the impact of residuals (that are unavoidably present in the 
considered maps) on the estimators for alignments of low multipoles. Uncertainties can 
be potentially introduced in each step of data analysis needed to generate the CMB map. 
For example, it is known that it is difficult to perform an accurate component separation 
close to the Galactic plane. Therefore we propose a new consistency test that can be 
capable of detecting the effects of the residuals on the estimators for the alignments at 
low multipoles. The basic idea of this check is the comparison between the value of the 
estimator of the map under analysis with the values of the estimator of the simulated maps 

^http://lambda.gsfc. nasa.gov/ 
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that are unlikely (or likely, depending on the context) at the same level of the considered 
map. The application of the proposed test to the not cleaned V band map, certainly not 
meaningful for a cosmological analysis, is performed in order to verify the capability of the 
method to identify spurious contributions in a known dirty case. 

The problem of minimizing potential residuals in a CMB map is crucial for cosmological 
analyses. In [43] it is proposed a way to limit the impact of foreground contamination 
present in the all sky maps adopted for the low I analysis of the alignments. This new 
method makes use of a power equalization filter [44] and it has been proposed for a better 
control of residual foregrounds and therefore for a potentially more robust cosmological 
analysis. However, even if this (or some other) component separation method could work 
perfectly, residuals can still be present because of uncertainties from potential systematics 
and from other stages in the data analysis. In [45] the a£„j (coefficients of the expansion 
over the basis of the Spherical Harmonics) of the CMB are obtained with a minimum 
variance method in case of incomplete sky coverage, non uniform noise and foreground 
contamination. Unfortunately, if the considered mask excludes more than about 10% of 
the sky then the error bars of the coefficients become too large. 

In this paper we adopt a drastic approach and mask the Galactic region (with masks of 
various size '^). Since the Multipolc Vectors expansion (that is the mathematical tool that 
we will use to define alignments [7, 8]) is doable only over the all sky (i.e. over the full set of 
the orthonormal Spherical Harmonics), we do not look for a new basis defined in the uncut 
region but we screen the information coming from the Galactic region, setting the value 
of the pixels falling into the area we do not want to consider, to a fiducial value (namely 
0). In Section 2 it is shown how to consistently perform the analysis and how to use the 
masking process in order to analyze the contamination present in the masked area (i.e. the 
proposed test) Technically, we study the redistribution (due to the masking process) of 
each bin of the probability distribution functions of the alignment estimators. The used 
technique can indeed only exclude that the considered map is clean in some region but 
it cannot state that the map is clean in that region since a possible residual could be in 
principle compatible even with random realizations. By construction, our method works 
for any possible mask suited to exclude regions that are possibly affected by various kinds 
of contamination. 

The paper is organized as follows: in Section 2 we present the adopted methodology 
and describe the performed simulations and the proposed test, in Section 3 we show the 
obtained results, and in Section 4 we draw our conclusions. 

2. Methodology 

2.1 Multipole Vectors 

The alignment of multipoles can be defined using a new representation of CMB anisotropy 

^In 2 cases out of 3, larger than 10% of the sphere. 

*As it will be clearer in Section 2, we do use the masking process as a tool to define a test for residuals 
and to infer conclusions about the all sky maps but we do not use the scientific analysis on the masked 
maps to infer conclusion about the all sky maps. 
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maps where the are replaced by vectors [7, 8]. In particular, each multipole order i is 
represented by £ unit vectors and one amplitude A 



(2.1) 



Note that the number of independent objects is the same in the l.h.s and r.h.s. of equation 
(2.1): 2£+l for aim equals 3£ (numbers of components of the vectors) +1 (given by A^^^) 
—i (because there are i constraints due to the normalization conditions of the vectors). 

Equation (2.1) can be understood starting from this observation [8]: if / is a solution 
of the Laplace equation 

VV = 0, (2.2) 

where = d'^ + dy + d'^ in Cartesian coordinates, then it is possible to build a new solution 
/' applying a directional derivative to / 

V^f = u-Vf = f', VY = 0, (2.3) 

with the gradient V = {dx,dy,dz)- This happens because the two operators and V^;: 
commute. Maxwell [46] repeated this observation £ times considering the 1/r potential as 
starting solution. Here r = {x,y,z) and r = Vf- f = ^/x^ + y'^ + z^. In this way, one 
obtains 

/^(x, y, z) = Vs, • • • Va, Vff, - . (2.4) 

Observe the simple pattern that emerges as we apply the directional derivatives one at a 
time: 



/o 
/i 
h 
h 



(-l)(^i-r) 



(3 • l)(tti • f)(tt2 • r) + r^(-ui ■ U2) 



(-5 • 3 • l)(ni • r){u2 ■ r)(^?3 • ^ + r'^{-) 



r7 



The (...) stands for a polynomial which we do not write explicitly, being not relevant to 
the current purposes. 

Moreover, writing ji in spherical coordinates once r is set to 1, one finds the following 
property 

VVK1,^,'^)=^(^+1)/^(1,^,0), (2-5) 
where is the angular Laplace operator defined as 



V2 = - 



sm y ~ ' sm^ ' 



(2.6) 



-4- 



In other words /^(l, 9, 4>) is eigenfunction of the angular part of the Laplace operator with 
eigenvalue given by + 1). This is nothing but the definition of spherical harmonics Y^^^ 
(e.g. see [47]). Therefore, for every £ we can write 

e 

0, 0) = J2 (^emYemiO, 4>) , (2.7) 
m=-e 

where the amplitude A^^^ has been inserted because of normalization purposes. Equa- 
tion (2.7) makes evident the association represented by equation (2.1). From equation (2.7) 
it is possible to write down the set of equations that has to be solved to pass from a^rn to 
multipole vectors. In order to see that this set is solvable we count the equations and the 
unknowns involved in this set. Prom equation (2.7) we have 2^+1 equations (one equation 
for each independent ^) plus £ equations from the normality conditions of the vectors 
(i.e. Ui ■ Ui = 1 where i runs from 1 to £). Therefore the total number of independent 
equations is 31+1. This is also the number of unknowns because we have 3 unknowns for 
each vector plus 1 given by the amplitude A^^'> . This shows that the set is solvable. 

One of the advantage of Multipole Vectors representation is that from these unit vectors 
one can easily construct scalar quantities that are invariant under rotation. Note that is 
not equally easy to obtain scalar quantities directly from the aem coefficients since they 
depend on the coordinate system. For a more detailed explanation of equation (2.1) and 
of the properties of that association sec for example references [7, 8, 13, 26]. 

Unfortunately, an explicit analytical expression for the association given in equation 
(2.1) is possible only for £ = 1. For £ ^ 1 numerical methods are needed ^. The Copi et 
al.'s algorithm (which use is acknowledged here) for constructing multipole vectors from a 
standard spherical harmonic decomposition is described in [7] and the implementation of 
it is public available ^. Other methods exist but, as far as we know, their implementation 
is not public available on a standard platform (see for example [8, 50] where the problem 
of finding i vectors is translated into the problem of finding the zeros of a polynomial of 
degree 2£). 

2.2 Alignment estimators 

We focus on the alignment quadrupole-octupole. Therefore we consider the following esti- 
mators widely used in literature (e.g. [8, 13, 51, 10, 7]): 

3 

523 = ^|g-o;|/3, (2.8) 

i=l 

3 

D23 = ^\q- di\/3. (2.9) 

i=l 

®In fact we would have 4£ + 1 equation because each i different from has a real and imaginary part. 
But considering that aem with m > are related to those with m < 0) through aj^ = {—l)"^ae-m we are 
left with 2£ + 1 equations (i.e. the expression must be real). 

^Indeed, for ^ = 2 it is possible to obtain the multipole vectors computing the eigenvectors of a symmetric 
and traceless tensor representing the quadrupole, see [48, 49]. 

'^http:/ / www.phys.cwru.edu /projects / mpvectors / 
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Here the symbol " stands for a vector with norm equal to 1. The "area vectors" are defined 
as 



q = q2i X q22 , 

01 = 632 X 633 , 

02 = 033 X 631 , 

03 = 031 X 632 , 



(2.10) 
(2.11) 
(2.12) 
(2.13) 



where q2j represent the two normalized multipole vectors (j = 1, 2) associated to the 
quadrupole, 63^ represent the three normalized multipole vectors (j = 1,2,3) associated to 
the octupole. 

Notice that all the estimators belong to the interval [0, 1] and contain absolute values 
in order to make them invariant under the reflection symmetry (see for example [13, 26]). 

2.3 Description of the performed simulations and test for residuals 

We have performed 3 x 10^ Gaussian random extractions of ACDM skies. We have masked 
each extraction ^ setting the pixels inside the mask to the fiducial value equal to 0. This 
allows us to screen the information present in the Galactic region and, at the same time, 
to have an all-sky map (even if artificial) such that the multipole vectors decomposition 
is still doable ^. We have largely exploited the extended temperature mask available at 
LAMBDA web site that has a sky coverage of 73% and other two masks whose percentage 
of sky coverage are around 83% and 93% For each realization (masked and not masked) 
we have computed the estimators defined in Subsection 2.2 passing from to multipole 
vectors through the Copi et al.'s routine. This allows us to compute the probability dis- 
tribution function (henceforth pdf ) of the considered estimators in the all-sky case and in 
the masked cases. In all the cases (masked and unmasked) we have computed the value 
of the estimators for four maps: the ILC map, the MCMC map, the V band map and the 
foreground reduced V band We do not use the masked pdf to directly extract scientific 
conclusions because this would not tell anything about the all sky map that is what we are 
interested in. What we do use is the masking process itself. Here it is a description of the 
test we propose. It consists of the following steps: 

• consider the random maps whose all-sky estimator value belong to the same bin of 
the pdf as the all-sky value of the estimator of the map under analysis; 

• mask these maps (both the set of random maps and the map under analysis that 
belong to the same bin) as described above; 

^Hereafter we will use random extraction, random realization and random map as synonyms. 

®In this paper wc will call masked maps, the maps that present the value for the pixels that fall into 
the mask, even if, strictly speaking, they are still all-sky. 

^"Wo use "sky coverage" to mean "observed sky" (that has not to be confused with masked area). 

^^A part the V band map that is taken into account as test case, all the other maps are only CMB maps 
obtained with different techniques. That's why we compare the same simulations with all these maps. 
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Figure 1: MoUweide projections of the considered maps. ILC map (upper left), MCMC map 
(upper right), foreground reduced V band map (lower left), V band map (lower right). The number 
of pixels in each map is Npix = 3145728 except for the MCMC map where Npi^ = 49152. A part 
the V band map that is taken into account as test case, all the other maps are only CMB maps 
obtained with different techniques. See also the text. 

• transform to the multipole vectors representation and compute the masked values of 
the estimator (both for the set of random maps and for the map under analysis that 
belong to the same bin); 

• build a pdf of these masked values obtained from the random extractions (what we 
call redistribution); 

• compare the masked value of the estimator for the map under analysis with the 
redistribution. 

In other words, we compute how each bin of interest of the all-sky pdf (for the considered 
estimator, D23 and S23 in this work) is redistributed when a mask is applied. This allows 
us to analyze how consistent are the four considered maps with the random realizations 
whose D23 or S23 values belong to the same bin of the all-sky distributions. More explicitly, 
if the estimator of the masked map under analysis does not behave consistently with the 
other random masked maps then we conclude that some residual was present in the all-sky 
map affecting the value of the all-sky estimator and making it belonging to a "wrong" 
bin. This consistency check is the test we propose in this paper to look for the effects of 
residuals, that are potentially present in the Galactic region of the considered WMAP Syr 
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maps, on the alignment estimators. Of course, in the case of the not cleaned V band map 
we are sure that a non-CMB signal is present at low Galactic latitudes. 

Note that the process adopted in our test introduces a sort of kernel representing the 
effect of the mask for each considered bin. The global pdf in the presence of a mask can 
be seen as the convolution of the all-sky pdf with this kernel. 

The resolution that we have considered to perform this analysis is represented by the 
HEALPix [52] parameter A^^jfie that we have set to 16. 

The impact of the noise in our test is addressed at the end of the next section. 

3. Results 

In this Section wc present and discuss the obtained results. We have chosen the following 
convention, in order to display our results on the pdf: red color for the all-sky pdf, cobalt 
blue color for the masked pdf, green color for the difference between masked and unmasked 
pdf. The green plots give the effect that has been introduced by the considered treatment 
of the mask. More precisely they show the difference between masked and unmasked 
distribution (masked minus unmasked). Moreover we use the blue color for the pdf of the 
redistribution for the large mask, the light blue color for the pdf of the redistribution for 
the medium sized mask and the white color for the pdf of the redistribution for the small 
mask. Vertical lines show the values of the estimators derived from the considered maps: 
black vertical line is for ILC map, blue vertical line for the MCMC map, green vertical 
line for the foreground reduced V band, yellow line for V band map. The V band map has 
been taken into account just for comparison with the foreground reduced V band and for 
method validation. 

3.1 Quadrupole-Octupole alignment estimators 

In Fig. 2 we show the all sky pdf for D23 with the corresponding values for the considered 
maps. In Fig. 3 we show the same estimator D23 for the three considered masks (left 
column of panels). We note that for all the considered maps, but not for the V band map 
(yellow line) where foreground subtraction is not applied, the D23 estimator assumes 
similar values when intermediate or large masks are applied, as the black, green, and blue 
lines tend to overlap (or at least to be closer) when the size of the Galactic mask increases. 
In the right column of panels of Fig. 3 there is the difference of the pdf 's between masked 
and all-sky case, that turns out to be stronger for larger masks. At the chosen binning 
of 0.01, such an effect appears weaker for the smaller mask It is interesting to notice 
that the area to the right of the black line (i.e. ILC map) below the pdf is larger for larger 

^^^httpi/Zhealpix.jpl.nasa.gov 

^^For the reader who is not famihar with HEALPix convention, Nside = 16 corresponds to maps of 3072 
pixels on the whole sky. 

^*The same is true also for S23, see the following and Fig. 6. 

^^Increasing the number of realizations (for example to 10®) or the size of the bin (for example to 0.05) 
would show a similar shape as the other masked cases. We do not report here these plots but we dedicate 
a short appendix to this aspect. 
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Figure 2: All-sky pdf for D23. Vertical lines represent the values for considered maps (of course 
in the all sky case). More precisely, black vertical line is for ILC map, blue vertical line for the 
MCMC map, green vertical line for the foreground reduced V band, yellow line for V band map. 
The panel presents the counts (y-axis) versus the statistic (x-axis). See also the text. 

masks. In particular, the probability to obtain a value smaller than the ILC-WMAP value 
is 98.36% for the all-sky case, 97.96% for smaller mask, 88.62% for the medium mask, 
87.37% for the extended mask. This would seem to indicate that excluding/screening 
the information coming from the Galactic region (as described in Section 2.3) makes the 
estimator no longer anomalous (for the map under analysis). Unfortunately this does not 
indicate uniquely that the anomaly that is present in the all-sky ILC-WMAP map is due 
to residuals that are present in the Galactic region. The fact that the masked ILC-WMAP 
map is no longer anomalous might also be an effect of the masking process itself (i.e. of the 
prior of setting to the masked pixels) To discriminate between the two possibilities 
we compute the redistribution of the 87*^ bin (that is the all-sky bin where the D23 value 
for the ILC-WMAP map falls) caused by the presence of zeros of the masking process. We 
analyze then if the masked ILC-WMAP D23 value is consistent with this redistribution, 
i.e. with a masked Gaussian random realization that exhibits a value for D23 that falls in 
same bin as the all-sky ILC-WMAP (this is the consistency check we propose). Of course, 
the same has been done for the other maps and for all the three considered masks. Fig. 4 
shows the "masking flux" of the bin 87 for the ILC-WMAP, of the bin 50 for the MCMC 
map, of the bin 96 for V band map and of the bin 26 for the foreground reduced V band 
map. Observing the fourth column of Fig. 4 we can see that the ILC-WMAP map is always 
consistent with a Gaussian random realization that shows the same all-sky D23 value. For 
what concerns the other maps. Fig. 4 shows that they are not consistent with a Gaussian 
random realization because first, second and third column of Fig. 4 present the D23 value 
for the other three maps (vertical lines) far from the peak of the distribution. We interpret 
this as evidence for the presence of residuals in the Galactic region in the all starting all-sky 
maps (except the aforementioned ILC-WMAP map). 

^^Therefore, as aforementioned in footnote 4 and in Subsection 2.3, we cannot extract scientific informa- 
tion about tlie all sky map from the masked analysis. 
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D23 D23 

Figure 3: Left column of panels: pdf of D23 for the masked case. Black vertical line is for ILC 
map, blue vertical line for the MCMC map, green vertical line for the foreground reduced V band, 
yellow line for V band map. Right column of panels: difference of the pdfs between masked and 
unmasked distributions. Percentage of sky coverage from top to bottom row of panels: 73%, 83% 
and 93%. All the panels present the counts (y-axis) versus the statistic (x-axis). See also the text. 

In Fig. 5 we show the all sky pdf for S23 with the corresponding values for the consid- 
ered maps. In Fig. 6 we show the same estimator S23 for the three considered masks. In 
particular, the probability to obtain a value smaller than the ILC-WMAP value is 98.90% 
for the all-sky case, 99.24% for smaller mask, 98.50% for the medium mask, 95.44% for 
the extended mask. As for the D23 estimator, this does not necessarily mean presence 
of residuals coming from the Galactic region. In Fig. 7 we show the redistribution for 
the bins where the S23 values of the considered maps stand. In analogy with what has 
been obtained for the redistribution of the D23 bins, Fig. 7 shows that the ILC-WMAP 
map is the only map that is consistent with the corresponding pdf of the redistribution 

^^In fact the V band for the larger mask might be considered consistent, but the inconsistency become 
evident as soon as it is considered a less sized mask. See also section 4. 
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Figure 4: Pdf of the redistribution of some specific bin of D23. In the first column we plot the 
redistribution due the masking process of the bin 26 for the foreground reduced V band map (green 
vertical line). In the second column we plot the redistribution due the masking process of the bin 
96 for the V band map (yellow vertical line). In the third column we plot the redistribution due 
the masking process of the bin 50 for the MCMC map (blue vertical line). In the fourth column 
we plot the redistribution due the masking process of the bin 86 for the ILC map (black line) . The 
first row is for the largest mask (73% of sky coverage), the second row for the middle size mask 
(83% of sky coverage) and the third row is for the smallest considered mask (93% of sky coverage). 
All the panels present the counts (y-axis) versus the statistic (x-axis). See also the text. 



Note that the pdfs of redistributions of Figs. 4 and 7 are more peaked as larger is the sky 
coverage. 

The redistributions shown in Figs. 4 and 7 have well defined shapes. The case of V 
band map and D23 estimator is the most anomalous (see Fig. 2), thus the corresponding bin 
is less populated (282 points). In spite of this, its redistributions are pretty stable. In order 
to show this, we take the case of the large mask and compute in Fig. 8 the redistribution 
for all the available points (left panel, as in Fig. 4) and for half of the available points 
(right panel). Since the shape appears to be stable for this case, the most critical one, we 
conclude that the shown results are robust. In general, when the number of realization are 
few, an increase of the number of realizations and/or of the bin size of the histogram of 
the redistribution might be needed. 

Regarding the comparison between our estimators for the V band map and the fore- 
ground reduced V band map, we note that the analysis of the pdf of the redistribution in 
both cases always indicates the presence of residuals in the Galactic region, except perhaps 
for the case of the S23 estimator applied to the V band map as suggested by the largest 
mask case (see Fig. 7). Note that the power of our method depends also on the chosen 
estimator and on the size of the mask. For this reason, given the estimator, we recommend 

^^This exception is not so surprising since the nature of our test (note also that the CMB temperature 
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Figure 5: The same as Fig. 2 but for S23. Vertical lines represent the values for considered maps 
(of course in the all sky case). More precisely, black vertical line is for ILC map, blue vertical line 
for the MCMC map, green vertical line for the foreground reduced V band, yellow line for V band 
map. The panel presents the counts (y-axis) versus the statistic (x-axis). See also the text. 

to apply our test exploiting various mask sizes. Regarding the power of the method, in the 
considered cases, we found that the width of the pdf of redistribution increases (although 
not too strongly) with the size of the mask. Note that for a map where we are aware 
of the presence of foregrounds, as the V band map, the test successfully detects spurious 
contamination. 

We have also studied the impact of the noise in the pdf of S23 and D23 in the masked 
(73% of sky coverage) and all-sky case. We expect only a (very) weak impact of the noise 
for our analysis as intuitively suggested by the fact that at Ngid^, = 16 the root mean 
square (RMS) of the signal (i.e. CMB) is ~ 50 and the RMS of the noise is ~ 2 /jK 
(for example considering the V band map of WMAP). In Fig. 9 we show the effect of the 
noise on the pdf confirming the weak impact. These plots have been obtained with 300000 
Gaussian random realizations. The CMB signal is extracted according to a ACDM model 
and the Gaussian random noise has been generated exploiting the non-uniform sensitivity 
of the V band map. 

4. Conclusions 

We have studied how the alignment between quadrupole and octupole is polluted on the 
Galactic region. 

We found that among the considered maps the WMAP-ILC is the most consistent 
with Gaussian random realizations. Prom the point of view of this test, we can say that 
WMAP-ILC is the most clean map (among the considered ones). This result has been found 
studying the redistribution (due to the masking process) of each bin of the probability 
distribution functions of the alignment estimators. Of course we cannot exclude that 
other tests might detect some residual contamination. The used technique can indeed only 

anisotropy map in tlie V band is clearly not so affected by diffuse Galactic foregrounds wiien a large fraction 
of the sky at low Galactic latitudes is excluded). 
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S23 S23 

Figure 6: The same as Fig. 3 but for S23. Left column of panels: pdf of S23 for the masked case. 
Black vertical line is for ILC map, blue vertical line for the MCMC map, green vertical line for 
the foreground reduced V band, yellow line for V band map. Right column of panels: difference of 
pdfs between masked and unmasked distributions. Percentage of sky coverage from top to bottom 
row of panels: 73%, 83% and 93%. All the panels present the counts (y-axis) versus the statistic 
(x-axis). See also the text. 

exclude that the considered map is clean in some region but it cannot state that the map is 
clean in that region since a possible residual could be in principle compatible with random 
realizations. 

Since the WMAP-ILC passes this test for three Galactic masks, the unlikely align- 
ment between Quadrupole and Octupole (probed by the two considered estimators) that 
is present in this map is confirmed. 

Although the results presented here formally apply to the considered products, our 
study points out on the relevance of analyzing low multipole alignments with the possibility 
of testing sky regions potentially affected by systematics (of instrumental or astrophysics 
origin) and show how it is feasible through Monte Carlo simulations. By construction, 
our method works for any possible mask adopted in the analysis of data from current and 
forthcoming CMB anisotropy experiments, such as those expected by the Planck satellite. 
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Figure 7: Pdf of the redistribution of some specific bin of S23. In the first column we plot the 
redistribution due the masking process of the bin 21 for the foreground reduced V band map (green 
vertical line) . In the same second column we plot the redistribution due the masking process of the 
bin 53 for the V band map (yellow vertical line) . In the third column we plot the redistribution due 
the masking process of the bin 48 for the MCMC map (blue vertical line). In the fourth column we 
plot the redistribution due the masking process of the bin 73 for the ILC map (black vertical line). 
The first row is for the largest mask (73% of the sky coverage) , the second row for the middle size 
mask (83% of the sky coverage) and the third row is for the smallest considered mask (93% of the 
sky coverage). All the panels present the counts (y-axis) versus the statistic (x-axis). See also the 
text. 
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A. Comparison of simulations with more realizations 

In Fig. 10 we focus on the smallest mask case. This case for the difference of pdfs seems 
not to have a clear shape (see lower-right panels of Fig. 3 and 6). Therefore we increased 
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Figure 8: Convergency of our results. Redistribution of D23 estimator, for bin 96 for the large 
mask case (73% of sky coverage). Left panel is exactly the same as in Fig. 4, while right panel uses 
half of the realizations belonging to the same bin. The yellow vertical line represents the value of 
the estimator for the V band map. See also the text. 



si 

TS 
ft 



rfi 100 



- all- sky 




IM 















Figure 9: Impact of the noise. Upper panels: difference of the pdf's for the D23 estimator. Lower 
panels: difference of the pdf's for the S23 estimator. Left panels: difference of the pdf in the case 
with noise minus without noise in the all-sky case. Right panels: difference of the pdf in the case 
with noise minus without noise in the masked case (extended mask, i.e. 73% of sky coverage). All 
the panels present the counts (y-axis) versus the statistic (x-axis). See also the text. 



the number of extractions in order to improve the statistics. In Fig. 10 we show the results 
for differences of the pdf's between masked and unmasked distribution of D23 and S23 
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Figure 10: Impact of the number of realizations on the considered estimators for the smaller mask 
case (93% of sky coverage). All the panels are referring to the smaller mask case (93% of sky 
coverage). Upper panels: difference of the pdfs between masked and full sky of D23 for the masked 
case. Lower panels: difference of the pdfs between masked and full sky of S23 for the masked 
case. Left panels: 1000000 extractions with binning of 0.01. Right panel: 1000000 extractions 
and binning of 0.05. See lower-right panels of Figs 3 and 6 for comparison with the corresponding 
pdf built with 300000 extractions. All the panels present the counts (y-axis) versus the statistic 
(x-axis). See also the text. 



with 1000000 extractions with binning of 0.01 and 0.05. This demonstrates that even the 
smaller mask case exhibits qualitatively the same behavior of the other cases (even if weaker 
quantitatively). Anyway, we suggest to use larger number of realizations when small masks 
are exploited. 
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